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Abstract 

Simple Temporal Networks (STNs) provide a powerful and general tool for representing con¬ 
junctions of max;imum delay constraints over ordered pairs of temporal variables. In this paper we 
introduce Hyper Temporal Networks (HyTNs), a strict generalization of STNs, to overcome the 
limitation of considering only conjunctions of constraints but maintaining a practical efficiency in 
the consistency check of the instances. In a Hyper Temporal Network a single temporal hyperarc 
constraint may be defined as a set of two or more maximum delay constraints which is satisfied 
when at least one of these delay constraints is satisfied. HyTNs are meant as a light generalization 
of STNs offering an interesting compromise. On one side, there exist practical pseudo-polynomial 
time algorithms for checking consistency and computing feasible schedules for HyTNs. On the 
other side, HyTNs offer a more powerful model accommodating natural constraints that cannot 
be expressed by STNs like “Trigger off exactly 5 min before (after) the occurrence of the first 
(last) event in a set.”, which are used to represent synchronization events in some process aware 
information systems/workflow models proposed in the literature. 

Keywords. Simple Temporal Networks Temporal Consistency Hypergraphs Pseudo-polynomial 
Time Algorithms Mean Payoff Games Disjunctive Temporal Problems Workflows 


1 Introduction 

In many areas of Artificial Intelligence (AI), inclnding planning, schednling and workflow manage¬ 
ment systems, the representation and management of qnantitative temporal aspects is of crncial im¬ 
portance [4,10,11,20,37,43]. Examples of possible qnantitative temporal aspects are: constraints on 
the earliest start time and latest end time of activities, constraints over the minimnm and maximnm 
temporal distance between activities, etc. 

In many cases these constraints can be represented as an instance of a Simple Temporal Network 
(STN) [18], a directed weighted graph where each node represents a time-point variable (timepoint), 
nsnally corresponding to the beginning or the end of an activity, and each arc specifies a binary 
constraint on the schednling times to be assigned to its endpoints. In [18], each arc is labeled with 

a closed interval of real valnes: for example, the labeled arc u v encodes the binary constraint 
X < V — u < y over its endpoints u and v. A more nniform and elementary representation of an STN 
is provided by its distance graph^ [18], a graph having the same set of nodes as the original one, bnt 

where each arc u ^ y fg replaced by two arcs, each labeled with a single real valne: arc u —^ v to 
express the constraint v — u < y, and arc v — A u to express the constraint u — v< —x, i.e., x < v — u. 

An STN is said to be consistent if it is possible to assign a real valne to each timepoint so that all 
temporal constraints are satisfied. The consistency property can be verified by searching for negative 
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^Distance graph is also called constraint graph by other authors |17]. Moreover, Bellman |3] was the first to describe 
the relation between shortest paths and difference constraints in a constraint graph. 
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cycles in the distance graph and it is well known that the consistency check and the determination of 
the earliest/latest value for each timepoint can be done in polynomial time [18]. 

However, STNs do not allow the expression of constraints like “trigger off an event exaetly 5 min 
after the oceurrence of the last of its predecessors”, which are a quite natural constraints to represent 
synchronization events in a process aware information system plan/workflow schema [25]. This is 
because in STNs, and in some of their natural extensions, (1) it is not possible to represent a single 
constraint involving more than two timepoints and (2) all constraints have to be satisfied in order to 
have the network consistent. On the contrary, the above constraint about a synchronization event can 
be represented as a set of distance constraints, each involving a different pair of timepoints, that is 
considered satisfied when at least one of set components is satisfied. In order to represent and analyze 
disjunctive constraints like the above one, it is then necessary to consider models like Disjunetive 
Temporal Problem (DTP) [44] where a constraint is a set of disjunctive difference constraints over the 
timepoints. The drawback of such model is that the consistency check problem is NP-complete [44]. 

1.1 Contribution 

In this article we propose to generalize STN to Hyper Temporal Network (HyTN), which allows also 
the expression of constraints like the above one regarding synchronization events, but where the 
consistency check is amenable of effective solution algorithms. 

Moreover, we show an interesting link between the consistency check of HyTNs and resolution in 
Mean Payoff Games (MPG), a family of perfect information games played on graphs by two oppo¬ 
nents [21], for which some pseudo-polynomial time algorithms for determining winning strategies are 
known [7,47]. 

A preliminary version of this article appeared in the proceedings of TIME symposium [15]. Here 
we extend the presentation as follows: (1) the definition of HyTN has been extended in order to 
allow the presence of two kinds of hyperarcs; (2) the motivating example section has been revised 
to show how the new kind of hyperarc can be used; (3) some further issues and pertinent properties 
about HyTN have been introduced and proved; (4) several proofs have been expanded and clarified; 
(5) the experimental analysis of the consistency check algorithm has been improved considering more 
recent algorithms [7] for finding winning strategies for MPGs. This has improved the performances 
dramatically. 

Organization The remainder of the article is organized as follows. In Section 2 we present a 
motivating example from the domain of the workflow-based process management to bring out HyTNs. 
Section 3 introduces some definitions and well-known results for STNs and introduces some definitions 
about hypergraphs. The generalization of STNs into HyTNs and the definition of consistency problem 
for HyTNs are presented in Section 4. In Section 5 we recall the main facts and results about Mean 
Payoff Games which are useful for the following sections. Section 6 presents the investigation into the 
link between the HyTN consistency problem and Mean Payoff Games deriving pseudo-polynomial time 
algorithms for checking the consistency of HyTNs and computing feasible schedules whenever they 
exist. Some empirical evaluations of the proposed algorithms are reported in Section 7. In Section 8, 
some related works are presented and discussed with respect to our approach. Section 9 summarizes 
the main facts brought to light in this article and presents a possible future development of the work 
we are currently carrying on. 


2 Motivating Examples 

In the introduction we have briefly recalled a kind of constraint that cannot be expressed within STNs. 
In this section, we describe in more detail two examples of temporal constraints that cannot be fully 
described in an STN in order to introduce and motivate the new expressive capability of our model. 
As a further motivation, at the end of the section we also spotlight how this new capability has been 
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Figure 1: A simple workflow schema excerpt with three parallel flows of execution. 


recently exploited to check the consistency of Conditional Simple Temporal Networks (CSTNs) [45] 
in a more efficient way. 

Let us consider an example in the domain of the workflow-based process management, a domain 
concerned with the coordination and control of business processes using information technology. A 
workflow is a representation of a business process as the coordinated execution of activities by human 
or automatic executors (agents). A Workflow management system (WfMS) is a software system that 
supports the automatic execution of workflows [25]. In a WfMS, the management of temporal aspects 
is a critical component and in the literature there are many proposals on how to extend a workflow 
in order to represent and manage temporal constraints of a business process [4,8,10-13,19,20,23]. 
In particular, in [4,8,13,19,20] authors show how to represent and manage some kinds of temporal 
constraints using specific algorithms, while in [10-12,23] authors show how it is possible to represent 
and manage a wider class of temporal constraints exploiting models like Time Petri Nets [35] or STNs/ 
STNUs [36]. 

In this paper we consider an excerpt of the conceptual temporal model proposed by Combi et 
al. [ 12 ], where the specification of a temporal workflow is given by a workflow schema, a directed graph 
(also called workflow graph) where nodes correspond to activities and arcs represent control flows that 
define activity dependencies on the order of execution. Both nodes and arcs may be associated to 
temporal ranges to specify temporal constraints. There are two different types of activity: tasks and 
connectors. Tasks represent elementary work units that will be executed by external agents. Each 
task is graphically represented by a box containing a name and a temporal range that specifies the 
allowed temporal span for its execution. Connectors represent internal activities executed by the 
WfMS to achieve a correct and coordinated execution of tasks. They are graphically represented by 
diamonds and, as with tasks, each of them has a temporal range that gives the temporal span allowed 
to the WfMS for executing it. Every arc has a temporal property that gives the allowed times that 
can be spent by the WfMS for possibly delaying the consideration of the next activity after the end 
of the previous one. There are different kinds of connector that allow one to modify a control flow. 
Split connectors are nodes with one incoming arc and two or more outgoing arcs: after the execution 
of the predecessor, (possibly) several successors have to be considered for the execution. The set of 
nodes that can start their execution is given by the kind of split connector. A split connector can be: 
Parallel, Alternative or Conditional. Join connectors are nodes with two or more incoming arcs and 
one outgoing arc only. The types of activities considered in [12] are a subset of the possible activities 
specified by the Workflow Management Coalition [1,25[. 

Fig. 1 shows a simple workflow schema where the Parallel connector +i splits the flow into three 
parallel flows of execution (one for the sequence of tasks Ti and T 2 , one for task T 3 , and one for task 
T 4 and T 5 ) that have to be joined (synchronized) by the AND join connector +2 before continuing 
the execution; all temporal ranges are in minutes. 

Let us consider the connector + 2 ; according to the recommendations from the Workflow Manage¬ 
ment Coalition (WfMC) [25] and the temporal specification from [12], the execution of this connector 
requires to wait all incoming flows and, after the last incoming flow, to wait a time according to the 
connector temporal range before following the outgoing arc. In other words, the incoming flows can 
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Figure 2: An STN representing temporal aspects of the workflow depicted in Fig. 1. The dotted region 
emphasizes, within the workflow excerpt, the connections to an AND join connector. 


arrive at different instants but only when the last one arrives, the connector has to be activated in 
order to continue with the execution. 

Combi et al. [12] proposed a method to translate workflow schemata to STNs/STNUs [36] in order 
to analyze and validate all temporal aspects in a rigorous way. As already noted in [9] and [10], such 
translation cannot specifically represent the behavior of an AND join connector, because the kind of 
constraints in an STN/STNU is limited. Therefore, in [10], the authors proposed an adjustment of the 
translation of an AND join connector introducing for each incoming arc of the connector a buffer node 
connected with some determined new arcs and assuming a reasonable but fixed execution algorithm 
for the STN. In more detail, let us consider Fig. 2 that depicts the representation of workflow of Fig. 1 
by means of an STN following partially the method described in [10] (without loss of generality, here 
we convert task constraints as STN arcs instead of STNU contingent ones because we are interested 
only in the AND join conversion). Each activity of the workflow is represented by two STN nodes, one 
to represent the begin timepoint, Bi, one for the end one, Ei, and temporal ranges in the workflow are 
represented by STN arc labels. Regarding the translation of the AND join node + 2 , nodes representing 
the task endings on parallel flows, Et 2 , Et^, and Et^, are connected to buffer nodes 61 , & 2 , and 63 
that allow the parallel flows to complete their execution following only their temporal constraints. 
Then, bi, & 2 , and 63 are connected to node B^^ (which represents the begin instant of the AND 
join connector) by temporal constraints [ 0 ,fi], [ 0 ,^ 2 ] and [ 0 , 13 ], where the values ti,t 2 , and are 
determined during the workflow-to-STN conversion as explained in [10]. 

Now, let us consider a possible execution scenario. If 61 , 621 and 63 occur all together at instant 20, 
then, following the proposed temporal semantics [ 12 ], the only possible instant value for B^^ must 
be 20 while the updated STN allows any value in the range [20,20 + minjfi,f 2 ,ta}]- In [10], the 
authors showed that the right value is always the lower bound of such extended range and, therefore, 
it is sufhcient to adopt an early execution strategy in order to choose the right value for timepoint 
B+ . 

In other words, the proposed translation has two drawbacks: (1) it requires some preliminary 
computations for determining ti,t 2 , and ts values, and (2) the resulting STN admits some solutions 
that are not admissible by the semantics of the AND join connector. 

To specifically represent the behavior of an AND join connector with respect to its predecessor time 
points without auxiliary conditions or analysis, it is necessary to introduce a new kind of constraint 
based on hyperarcs, as shown in Fig. 3. In the figure the multi-tail hyperarc A consists of three dashed 
arcs—called components —and replaces the arcs from bi {i = 1,2,3) to B^^ of Fig. 2. We say that 
a multi-tail hyperarc is satisfied if at least one of its components is satisfied. In Fig. 3 dashed arcs 
define the hyperarc A that is satisfied if B^^ is 0 distant from at least one time point among 61 , & 2 , 
and & 3 . Since B^^ is constrained to occur at the same instant or after each time point 61 , & 2 , and 63 
by the arcs between B^^ and bi, i = 1,2,3, the result is that to satisfy A it is necessary that B^^ 
occurs at the same instant of the last time point among 61 , 62 , and 63 , as required originally. In more 
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Figure 3: An augmented STN, that we call HyTN, where dashed arcs represent components of hyperarcs, 
a new kind of constraint. This HyTN improves the representation of the STN of Fig. 2. To emphasize the 
changes, here we have summed all arcs outside the dotted region. 



(a) An excerpt of a workflow schema 
containing a structured discriminator 
connector, "tu. 



(b) A possible representation of 
temporal aspects of a discrimi¬ 
nator connectors by means of a 
multi-head hyperarc. 


Figure 4: 


A structured discriminator connector and a possible representation of its temporal aspects. 


general, a multi-tail hyperarc is defined as a set of distance constraints (components) between some 
time points and a common end point. 

The use of hyperarcs allows also the representation of temporal aspect of other advanced connectors 
as, for example, the Structured Discriminator [1]. The Structured Discriminator connector provides 
a means of merging two or more distinct flows in a workflow instance into a single subsequent. In 
particular, it triggers the subsequent flow as soon as the the first incoming flow arrives. The arrival 
of other incoming flows thereafter have no effect on the subsequent flow. As such, the Structured 
Discriminator provides a mechanism for progressing the execution of a process once the first of a 
series of concurrent tasks has completed and according to the connector temporal range. Fig. 4a 
depicts an excerpt of a workflow schema containing a structured discriminator connector, ds, that 
joins three parallel flows. 

At the best of our knowledge, currently there are no proposals for the representation of temporal 
constraints of a discriminator connector in any temporal workflow model or process-aware information 
system [32]. Even exploiting the methodology proposed in [10], it is easy to verify that it is not possible 
to represent such connector as an STN because in a consistent STN all constraints have to be satisfied 
while here it is necessary to allow the possibility that only one constraint of a set has to be satisfied in 
order to specifically represent a discriminator connector. A possible way for specifically managing a 
discriminator connector consists in following the approach suggested by [10] for representing activities 
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and considering a variant of the multi-tail hyperarc, called multi-head hyperarc, for representing its 
temporal constraint, as depicted in Fig. 4b. In the figure there is a multi-head hyperarc A that connects 
the node representing the beginning instant of the discriminator activity to all nodes representing the 
end instant of the activities that precede the considered discriminator and are directly connected to 
it. In general a multi-head hyperarc is defined as a set of distance constraints (components) between 
one time point and some end points. We say that a multi-head hyperarc is satisfied if at least one of 
its components is satisfied. In Fig. 4b, dashed arcs define the hyperarc A that is satisfied if B_^ is 
at least 2 distant from Et^ or 1 distant from Et^ or 5 distant from Et^- It is sufficient that one of 
such previous nodes is executed and that the delay represented in the corresponding connecting arc 
is passed to execute Bj_, as required by the structured discriminator connector semantics. 

HyTNs are not only suitable for better representing temporal constraints originating from temporal 
workflow, but also for better representing more general temporal constraint networks like Conditional 
Simple Temporal Network [45]. 

A Conditional Simple Temporal Network (CSTN) is an enriched graph for representing and reason¬ 
ing about temporal constraints in domains where some constraints may apply only in certain condition 
settings (scenarios). Each condition in a CSTN is represented by a propositional letter whose truth 
value is observed in real time as the outcome of the execution of an observation time-point. An execu¬ 
tion strategy for a CSTN has to determine an execution time for each time-point guaranteeing that all 
temporal constrains that are significant in the resulting scenario are satisfied. An execution strategy 
can be dynamic in that its execution decisions can react to the information obtained from such obser¬ 
vations. The Conditional Simple Temporal Problem (CSTP) consists in determining whether a given 
CSTN admits a dynamic execution strategy for any possible combination of propositional outcomes 
happens to be observed over time. If such a strategy exists, the CSTN is said to be dynamically 
consistent (DC). 

Tsamardinos et al. [45] solved the CSTP by first encoding it as a meta-level Disjunctive Temporal 
Problem (DTP), then feeding it to an off-the-shelf DTP solver. Although of theoretical interest, this 
approach is not practical because the CSTP-to-DTP encoding has exponential size, and the DTP 
solver itself runs in exponential time. To our knowledge, this approach has never been empirically 
evaluated [26]. 

In [16], Comin and Rizzi proposed a novel representation of CSTNs in terms of HyTNs allowing the 
determination of the first singly exponential-time algorithm for checking the dynamic consistency of 
Conditional Simple Temporal Networks. More precisely, a CSTN instance is represented as a suitable 
HyTN where each possible scenario is represented and connected to other scenarios in an appropriate 
way and, then, such HyTN instance is solved in pseudo-polynomial time by the algorithms analyzed 
in the present paper. 

In summary, HyTNs allow the representation of temporal constraints that are more general of those 
represented in STNs [18], because they allow disjunctions involving more than two time points, but 
less general than those represented in DTPs [44] because all disjunctions related to a multi-head(tail) 
hyperarcs have to contain a common variable. Such kind of STN generalization not only allows the 
compact representation of some common temporal constraints in the domains like the workflow-based 
process management but also allows the determination of new interesting algorithm for checking 
dynamic-consistency in richer models like CSTN. 


3 Background and Notation 

In this section, we introduce some definitions, notations and well-know results about graphs and 
conservative graphs; moreover, we recall the relation between the consistency property of STNs and 
the conservative property of weighted graphs. 

We consider graphs that are directed and weighted on the arcs. Thus, if G = (V, A) is a graph, 
then every arc a € A is a triple {ta,ha,Wa) where ta € V is the tail of a, ha € V is the head of a, 
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and Wa G R is the weight of a. Moreover, since we use graphs to represent distance constraints, they 
do not need to have either loops (unary constraints are meaningless) or parallel arcs (two parallel 
constraints represent two different distance constraints between the same pair of node: only the most 
restrictive is meaningful). We also use the notations h{a) for ha, t{a) for ta, and w{a) or w{ta, ha) for 
Wa, when it helps. The order and size of a graph G = {V,A) are denoted by n = \V\ and m = |A|, 
respectively. The size is a good measure for the encoding length of G. 

A eycle of G is a set of arcs GCA cyclically sequenced as oq, oi,..., ai-i so that h{ai) = t{aj) 
if and only if J = (i + 1) mod it is called a negative cycle if w{C) < 0, where w(C) stands for 
Z^eeC We- A graph is called conservative when it contains no negative cycle. 

A potential is a function p : 1^ i—R,. The reduced weight of an arc a = (u, v, Wa) with respect to a 
potential p is defined as wg = Wa — Pv + Pu- A potential p oi G = {V, A) is called feasible if > 0 
for every a € A. Notice that, for any cycle G, wP{G) = w{G). Therefore, the existence of a feasible 
potential implies that the graph is conservative as w{G) = w^{G) > 0 for every cycle G. 

The Bellman-Ford algorithm [17] can be used to produce in 0{nm) time: 

• either a proof that G is conservative in the form of a feasible potential function; 

• or a proof that G is not conservative in the form of a negative cycle G in G. 

When the graph is conservative, the smallest weight of a walk between two nodes is well defined, 
and, fixed a root node r in G, the potentials returned by the Bellman-Ford algorithm are, for each 
node V, the smallest weight of a walk from r to v. Moreover, if all the arc weights are integral, then 
these potentials are integral as well. Therefore, the Bellman-Ford algorithm provides a proof to the 
following theorem. 

Theorem 1 ([3,17,22]). A graph admits a feasible potential if and only if it is conservative. Moreover, 
when all arc weights are integral, the feasible potential is an integral function. 

An STN can be viewed as a weighted graph whose nodes are timepoints that must be placed on 
the real line and whose arcs express mutual constraints on the allocations of their end points. An 
STN G = (V, A) is called consistent if it admits a feasible scheduling, i.e., a scheduling s : i—>■ R 

such that 

s(v) < s{u) +w(u,v) V arc {u,v) of G. 

Corollary 1 ([3,17,18]). An STN G is consistent if and only if G is conservative. 

Proof. A feasible scheduling is just a feasible potential. Therefore, this corollary is just a restatement 
of Theorem 1. □ 

In this paper, we also deal with directed weighted hypergraphs. 

Definition 1 (Hypergraph). A hypergraph is a pair (V, A), where V is the set of nodes, and A is 
the set of hyperarcs. Each hyperarc A G A is either a multi-head or a multi-tail hyperarc. 

A multi-head hyperarc A = {tA, Ha,wa) has a distinguished node Ia, called the tail of A, and a 
nonempty set Ha C E \ {t/i} containing the heads of A; to each head v G Ha is associated a weight 
wa{v) G R. Fig. 5a depicts a possible representation of a multi-head hyperarc: the tail is connected 
to each head by a dashed arc labeled by the name of the hyperarc and the weight associated to the 
considered head. 

A multi-tail hyperarc A = {TA,hA,WA) has a distinguished node Ha, called the head of A, and 
a nonempty set Ta C E \ {hyi} containing the tails of A; to each tail v G Ta is associated a weight 
wa{v) G R. Fig. 5b depicts a possible representation of a multi-tail hyperarc: the head is connected 
to each tail by a dotted arc labeled by the name of the hyperarc and the weight associated to the 
considered tail. 

The cardinality of a hyperarc A G A is given by |A| = \Ha U {tyi}| if A is multi-head, and 
|A| = |r^ U {SaW if A is multi-tail; if |A| = 2, then A = {u,v,w) is a standard arc. The order and 
size of a hypergraph (V, A) are denoted by n = jEj and m = X]y4e.4 1^1’ respectively. 
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(a) Multi-Head Hyperarc 
A = {tA,HA,WA)- 



(b) Multi-Tail Hyperarc A = 
{TA,hA,WA)- 


Figure 5: A graphical representation of the two kinds of hyperarcs. 


4 Hyper Temporal Networks and Consistency Property 

We introduce now Hyper Temporal Networks (HyTNs), a strict generalization of STNs to partially 
overcome the limitation of allowing only conjunctions of constraints. Compared to STN distance 
graphs, which they naturally extend, HyTNs allow a greater flexibility in the definition of temporal 
constraints. 

A HyTN is a directed weighted hypergraph TL = {V,A) where a node represents a time point 
variable (timepoint), and a multi-head/multi-tail hyperarc represents a set of temporal distance con¬ 
straints between the tail/head and the heads/tails, respectively. 

For example, the multi-tail hyperarc A = (Ta, B^^,wa) in Fig. 3, where Ta = {&i,& 2 ,& 3 } and 
WA{bi) = 0 for f = 1, 2,3, stands for the set of distance constraints — bi < 0 \ i = 1,2, 3}. 

In general, we say that a hyperarc is satisfied when at least one of its distance constraints is 
satisfied. Then, we say that a HyTN is consistent when it is possible to assign a value to each 
time-point variable so that all of its hyperarcs are satisfied. 

More formally, in the HyTN framework the consistency problem is defined as the following decision 
problem. 

Definition 2 (General-HyTN-Consistency). Given a HyTN H = {V,A), decide whether there 
exists a scheduling s : H —>■ R such that, for every hyperarc A G A, the following holds: 

• if A = (t, h, w) is a standard arc, then 


s{h) — s{t) < w; 


• if A = {tA, Ha,wa) is a multi-head hyperarc, then 

s{tA) > min {s{v) - wa{v)}; 
vGHa 


• if A = {Ta, hA,WA) is a multi-tail hyperarc, then 

s{hA) < max{s(?;) (?;)}. 

vGTa 


Any such scheduling s : H —)■ R is called feasible. A HyTN that admits at least one feasible 
scheduling is called consistent. 

Comparing the consistency of HyTNs with the consistency of STNs, the most important aspect of 
novelty is that, while in a distance graph of a STN each arc represents a distance constraint and all 
such constraints have to be satisfied by a feasible schedule, in a HyTN each hyperarc represents a set 
of one or more distance constraints and a feasible scheduling has to satisfy at least one such distance 
constraints for each hyperarc. 
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Let us show some interesting properties about the consistency problem for HyTNs. 

The first interesting property is that any integral-weighted HyTN admits an integral feasible 
schedule when it is consistent, as proved in the following lemma. 

Lemma 1. Let 'H = (V, A) be an integral-weighted and consistent HyTN. Then % admits an integral 
feasible scheduling s : L {~T, —T -|- 1,..., T — 1, T}, where T = vCiV 

Proof. Since PL is consistent, then there exists a feasible scheduling s : L —?■ R. The idea in this 
proof is to project the HyTN PL over a conservative graph and then, in that setting, to exploit the 
integrality properties of potentials as stated in Theorem 1. However, this projection is asked to resolve 
the non-determinism contained in the disjunctive nature of the hyperarcs; in order to sort out such 
non-determinism, the projection is built considering the given feasible scheduling s as follows. 

For each hyperarc A G A, a weighted directed arc ca is defined as follows: 

• if H = {u,v,w) is a standard arc, then ca — {u,v,w). Note that s{v) < s{u) -f w follows by the 
feasibility of s; 

• ii A = {La, Ha, wa) is a multi-head hyperarc, then 

= {tA,VA,'WA{v)) where va = arg min {s(u) - wa{v)}. 

vGHa 


Here, s(u^) < s{tA) + WAiv) follows by the feasibility of s; 

• ii A = {Ta, La, wa) is a multi-tail hyperarc, then 

Ca — {vA,hA,WA{v)) where va = argmax{s(u) -|-w^(u)}. 

v€Ta 

Here, s(/i^) < s{va) + wa{v) follows by the feasibility of s. 

Now, a weighted directed graph G = (V, E) with E = {ca \ A G is defined. G is integral-weighted 
and conservative graph since it admits s as a potential function. Therefore, G admits an integral 
potential function s : H —>■ {—T,—T -|- 1,..., T — 1, T}. Indeed, such a function s is obtained by 
applying the Bellman-Ford algorithm on G. To conclude, we observe that s is also an integral feasible 
scheduling for PL. □ 

The following theorem states that General-HyTN-Consistency is NP-complete. 

Theorem 2. General-HyTN-Consistency is an NP-complete problem even if input instances 
PL = (V, A) are restricted to satisfy WAi') G {—1, 0,1} and \Ha\, \Ta\ < 2 for every A G A. 

Proof, li PL = {V, A) is integral-weighted and consistent, then it admits an integral feasible schedul¬ 
ing s : H —>■ {—r, ...,r} by Lemma 1. Moreover, any such feasible scheduling can be verified in 
polynomial time with respect to the size of the input; hence, General-HyTN-Consistency is in 

NP. 

To show that the problem is NP-hard, we describe a reduction from 3-SAT. 

Let us consider a boolean 3-CNF formula with n > 1 variables and m > 1 clauses: 

m 

ip{xi, . . . ,Xn) = /\{0H\/ PiV Ji) 

i=l 

where Ci = (ai V PiV ji) is the z-th clause of tp and each G {xj,Xj | 1 < j < n} is either a 

positive or a negative literal. 

We associate to p a HyTN PL^p = {V,A), where each boolean variable Xi occurring in p is rep¬ 
resented by two nodes, Xi and Xi. V also contains node z that represents the reference initial node 
for the HyTN PLp, i.e., the first node that has to be executed. For each pair Xi and Xi, PLp contains 
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[ 0 ] [ 0 ] 


(a) Gadget for a 3-SAT variable Xi. (b) Gadget for a 3-SAT clause Cj = (oj V 

pj V ) where each aj, pj , 7 ^ is a positive or 
negative literal. 

Figure 6: Gadgets used in the reduction from 3-SAT to General-HyTN-Consistency. 

a pair of hyperarc constraints as depicted in Fig. 6a: one with multi-head {xi,Xi} and tail in z and 
the other multi-tail {xi,Xi} and head in z. If is consistent, the pair of hyperarcs associated to 
X, ->x assures that admits a feasible scheduling s such that s{xi) and s{xi) are coherently set with 
values in {0,1} (see Lemma 1). In this way, s is forced to encode a truth assignment on the Xi’s. The 
HyTN TLip contains also a node Cj for each clause Cj of (p\ each node Cj is connected by a multi-tail 
hyperarc with head in Cj and tails over the literals occurring in Cj and by two standard and opposite 
arcs with node z as displayed in Fig. 6b. Such setting of arcs assures that if TL^p admits a feasible 
scheduling s, then s assigns value 1 at least to one of the node representing the literals connected with 
the hyperarc. 

More formally, T-L^, = {V, A) is defined as follows: 

• V = {z} U {xi I 1 < i < n} U {xi I 1 < f < n} U {Cj | 1 < j < m}; 

• A = Ur=i Vari U Cla^ , where: 

- Var* = ^{z,XiA),ixi,z,0),{z,XiA),(xi,z,0), 

{{Xi,Xi}, z,[w{xi),w{xi)] = [-1,-1]), 

(z,{xi,Xi}, [w(xi), r(;(xi)] = [0,0])|. 

This defines the variable gadget for Xi as depicted in Fig. 6a; 

- Claj = |(z,Cj, l),(Cj,z,-1), 
i{o:j,l3j,lj},Cj,[wiaj),w{l3j),w(jj)] = [0,0,0])|. 

This defines the clause gadget for clause Cj = (at V V ji) as depicted in Fig. 6b. 

Notice that |y| = 1 -I- 2n -I- m = 0{m + n) and rriA = 8 n + 5m = 0{m + n); therefore, the 
transformation is linearly bounded in time and space. 

We next show that (f is satisfiable if and only if T-Cp is consistent. 

Any truth assignment v : {xi,...,x„} —> {true, false} satisfying (p can be translated into a 
feasible scheduling s : F —s- Z of Hp as follows. For node z, let s(z) = 0, and let s{Cj) = 1 for each 
j = 1,..., m; then, for each i = 1,..., n, let s(xi) = 1 and s(Ti) = 0 if the truth value of Xi, ly(xi), is 
true, otherwise let s(xi) = 0 and s(xi) = 1. It is easy to verify that, using this scheduling s, all the 
constraints comprising each single gadget are satisfied and, therefore, the network is consistent. 
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Vice versa, assume that T-iip is consistent. Then, it admits an integral feasible scheduling s by 
Lemma 1. After the translation s{v) = s{v) — s(z), we can assume that s(z) = 0 . Hence, s(Cj) = 1 
for each j = 1,...,to, as enforced by the two standard arcs incident at Cj in the clause gadget, 
and {s(a;i), s(xi)} = { 0 , 1 } for each i = l,...,n, as enforced by the constraints comprising the 
variable gadgets. Therefore, the feasible scheduling s can be translated into a truth assignment 
ly : {a;i,..., Xn} {true, false} defined by iy{xi) = true if s{xi) = 1 (and s(Ti) = 0); iy{xi) = false 
if s{xi) = 0 (and s{xi) = 1 ) for every i = 1 , ..., n. 

To conclude, we observe that any hyperarc A € A of T-Lip has weights wa{-) G {—1, 0,1} and size 
|A| < 3. Since any hyperarc with three heads (tails) can be replaced by two hyperarcs each having at 
most two heads (tails), the consistency problem remains NP-Complete even if wa{-) G { — 1,0,1} and 
|A| < 2 for every A £ A. □ 

Theorem 2 motivates the study of consistency problems on hypergraphs having either only multi¬ 
head or only multi-tail hyperarcs. In the former case, the consistency problem is called Head-HyTN- 
CONSISTENCY, while in the latter it is called Tail-HyTN-Consistency. In the following theorem 
we observe that the two problems are inter-reducible, i.e., we can solve consistency for any one of the 
two models in /(m, n, W) time whenever we have a /(to, n, W) time procedure for solving consistency 
for the other one. 

Theorem 3. Head-HyTN-Consistency and Tail-HyTN-Consistency are inter-reducible by 
means of log-space, linear-time, local-replacement reductions. 

Proof. We show the reduction from multi-tail to multi-head hypergraphs; the other direction is sym¬ 
metric. Informally, what we will do is to reverse all the arcs (so that what was multi-tail becomes 
multi-head), and, contextually, we invert the time-axis (to account for the inversion of the direction 
of all arcs). 

Let Td = {V, A) be a multi-tail hypergraph, we associate to a multi-head hypergraph T-L' = (V, A') 
by reversing all multi-tail hyperarcs. Formally, we define 

A' = {(n, S, w) I {S, V, w) G A}. 

For example, a standard arc {t, h,w) G A is transformed into a reversed standard arc {h, t, w) in A! 
while a hyperarc with two weighted tails ti and t 2 becomes a hyperarc having ti and t 2 as its two 
weighted heads. 

Now, % is consistent if and only if T-L' is consistent. To prove it, we note that each scheduling s 
for % can be associated, with a flip of the time direction, to the scheduling s' = —s. Then, it holds 
that s is feasible for % if and only if s' is feasible for %'. Indeed, s satisfies the constraint represented 
by an hyperarc A = (Ta, Iia, wa) G A, namely 

s{hA) < max{s(n) -I- wa(w)}, 
vGTa 


or, equivalently 

if and only if s' (that is, 
ihA,TA,WA), namely 


-s(/ia) > min{-s(n) -^^(n)}, 
vGTa 

—s) satisfies the constraint represented by the reversed hyperarc A! = 
s'{H a) > min{s'('(;) —WA'iv)}- 


In the remainder of this work we shall adopt the multi-head hypergraph as our reference model. 
Hence, when considering hypergraphs and HyTNs, we will be implicitly referring to multi-head hy¬ 
perarcs. Notably, we consider the following specialized notion of consistency for HyTNs. 
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Definition 3 (HyTN-Consistency). Given a (multi-head) HyTN % = {V,A), decide whether there 
exists a scheduling s : G —>■ R such that: 

s(t^) > min {s(r)) — wa(i’)} 'iA&A. (1) 

v&Ha 

Remark 1. Notice that this notion of consistency for HyTNs is a strict generalization of STN one. 
In general, the feasible schedules of an STN are the solutions of a linear system and, therefore, they 
form a convex polytope. Since an STN may be viewed as a HyTN, the space of feasible schedules 
of an STN can always be described as the space of feasible schedules of a HyTN. The converse is 
not true because feasible schedules for a HyTN do not form a convex polytope. Let us consider, for 
example, a HyTN of just three nodes xi, X 2 , X 3 and one single hyperarc with heads {xi,X 2 } and tail 
X 3 expressing the constraint X 3 > min{a;i, 0 : 2 }; (0,2,2) and (—2,0,2) are both admissible schedules, 
but (1,1,0) = i(0, 2, 2) — i(—2, 0,2) is not an admissible schedule. In conclusion, the STN model is 
a special case of the Linear Programming paradigm, whereas the HyTN model is not. 

In the rest of this section, we extend the characterization of STN consistency recalled in Section 3 
to HyTNs. 

Definition 4 (Reduced Slack Value w\{v)). With reference to a potential p : V — )■ R, we define, for 
every arc A G A and every v G Ha, the reduced slack value w^{v) as wa(v) + p{tA) — p{v) and the 
reduced slack w\ as 

= max{rt;(^(u) | v G Ha}- 

A potential p is said to be feasible if and only if > 0 for every A G A. 

Again, as it was the case for STNs, a mapping f : V —)■ R is a feasible potential if and only if 
it is a feasible schedule. In order to better characterize feasible schedules, we introduce a notion of 
negative cycle. 

Definition 5 (Negative Cycle). Given a multi-head HyTN H = {V,A), a cycle is a pair {S,C) with 
S' C V and CCA such that: 

1. S = [J^^^{Ha U {Ia}) and S 0; 

2. \/v G S there exists an unique A G C such that Ia = v. 

Moreover, we let a(v) denote the unique arc A G C with Ia = v as required in previous item 2. Every 
infinite path in a cycle {S,C) contains, at least, one finite cyclic sequence Uj, ..., Uj+p, where 
Uj+p = Vi is the only repeated node in the sequence. A cycle (S, C) is negative if for any finite cyclic 
sequence vi,V 2 , ■ ■ ■ ,Vp, it holds that 


p-i 

< 0 . 

t=i 

There are two results about negative cycles as stated in the following lemmas. 

Lemma 2. A HyTN with a negative cycle admits no feasible schedule. 

Proof. By contraposition. Let 7^ be a consistent HyTN and let p be a feasible potential for H. Also, 
let (S, C) be any cycle of 77; we will show that (S, C) is not negative. For every A G C, let hA be the 
head of A with maximum reduced slack value: 

hA = arg max 
vSHa 
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Let us consider the infinite path in {S,C) built choosing, at each node Vt, ha(vt) the following 
node. As already seen, such a path contains at least one finite cyclic sequence Vh,Vh+i, ■. ■ ,Vk with 
Vk = Vh- The sum of weights of the finite cyclic sequence is given by 

fc-i fc-i 

t—h t—h 

for every potential p; since p is feasible, all terms of the last sum are non-negative. It follows that 
(S', C) is not negative. □ 

At first sight, it may appear that checking whether (S, C) is a negative cycle might take exponential 
time since one should check a possibly exponential number of cyclic sequences. The next lemma shows 
instead that it is possible to check the presence of negative cycle in polynomial time. 

Lemma 3. Let {S,C) be a cycle in a HyTN. Then checking whether {S,C) is a negative cycle can be 
done in polynomial time. 

Proof. Consider the weighted graph G = (S,UtgsAt), where each hyperarc At = {{t,v,—Wa(t){v)) \ 
V G iLa(t))}- 

Checking whether (S, C) is a negative cycle amounts to check whether all cycles in G have strictly 
positive weight. To do this, firstly, a potential function tt for G is determined by Bellman-Ford 
algorithm. If the algorithm returns a negative cycle instead of tt, then there is a negative cycle in 
{S,C) and the check ends. 

Otherwise, since w{G) = w^{G) > 0 for every cycle G of G, it is necessary to verify that no 
cycle in G has w'^{G) = 0. This control can be done by verifying the acyclicity of the subgraph of G 
comprising only arcs a of G with w'"{a) = 0. The check that a graph is acyclic can be done in linear 
time by a depth first visit [17]. □ 

A hypergraph "H is called conservative when it contains no negative cycle. In the next sections we 
will provide a pseudo-polynomial time algorithm that always returns either a feasible scheduling or a 
negative cycle, thus extending the validity of the classical good-characterization of STN consistency to 
general HyTN consistency. Here, we anticipate the statement of the main result in order to complete 
this general introduction of HyTNs. 

Theorem 4. A HyTN H is consistent if and only if it is conservative. Moreover, when all weights 
are integral, then H. admits an integral scheduling if and only if it is conservative. 

Proof. If H. is consistent, then it is conservative by Lemma 2. If 77 is not consistent, then there is a 
negative cycle as shown in Theorem 8-(3). The existence of an integral scheduling when all weights 
are integral is guaranteed by Lemma 1. □ 


5 Mean Payoff Games 

In this section, we propose an introduction to Mean Payoff Games (MPGs) tailored to the needs of 
the present work. MPGs represent a well-studied model for representing some kinds of two-player 
dynamics and we will show in Section 6 that there is a substantial equivalence between the MPG and 
the HyTN model, which will allow us to exploit some important algorithmic and structural results. 

An MPG is a weighted directed graph G = (Vq U ki, A) whose node set V is partitioned into two 
disjoint sets Vq and Vi, where, for p = 0,1, the nodes in Vp are those under control of Player p. Even 
with these graphs we have no loops and no parallel arcs. It is also assumed that every node has at 
least one outgoing arc. Notice that, in general, (Vq) Li) does not need to be a bipartiton of G, i.e., E 
may contain arcs with both endpoints in Vq, or with both endpoints in Vi. 

Each play starts with a pebble placed at some node uq G Lq U Vi and consists in a sequence of 
moves. Move t begins with the pebble placed in node vt-i and is played by the Player p such that 
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Vt-i G Vp'. the player chooses any arc e G E with tail te = Vt-i and moves the pebble along e; at the 
end of the move the pebble is in node Vt = he- The game ends as soon as Vt = Vt' for some t > t', i.e., 
when the pebble comes back to an already visited node Vt>. At this point, the pebble has traversed a 
cyclic sequence of arcs et/+i ,. ■. ,et and Player 0 “pays” to Player 1 the average weight of the visited 
cycle: 

i=t'+i 

If this amount is negative, then Player 0 wins the game, otherwise the winner is Player 1. 

A strategy for Player p is a mapping that, given all the previous visited nodes and the current 
node, returns which node has to be visited in the next move; a strategy is said to be positional (or 
memoryless) if it depends only on the current position vt and does not take into account all the 
previous history. If s € Vq U Vi and Player p has a strategy leading him to win any possible play 
starting at vq = s, then we say that s is a winning start position for Player p. We denote by Wp 
the set of winning start positions for Player p. A winning strategy for Player p leads Player p to win 
every play started from any node in Wp. Since these finite games are zero-sum, i.e., what won by a 
player is what lost by the other one, then they admit a game value v: for each start position s G P 
of the game, there exists a G R such that Player 0 has a strategy ensuring payoff at most Vg, while 
Player 1 has a strategy ensuring payoff at least Vg. 

It is worthwhile to consider an infinite variant of the model, in which the game does not stop, and 
continues for an infinite number of steps. In this model. Player 1 wants to maximize the limit inferior 
of the average weight: 

lim inf — 

n—¥oo Ti 

Symmetrically, Player 0 wants to minimize the limit superior of the same average weight: 

1 " 

limsup-y^ w{vt-i,vt) 

n—^cxD Tb 

In their Determinacy Theorem, Ehrenfeucht and Mycielski [21] proved that any infinite game 
admits a value and that this value equals the one of the finite counterpart game on every start 
position, i.e., = Vg for every s G VoUPi. Moreover, they proved the existence of positional strategies 

which are optimal for both variants of the model: when Player p limits himself to an optimal strategy 
TTp, i.e., when, for every v G Vp, he disregards all arcs with tail in v except the one with head in TTp{v), 
then he will secure himself the optimal payoff ly in every play, finite or infinite, however the adversary 
plays. The graph obtained from G by dropping all arcs with tail in Vp not prescribed by iTp is 
called the projection of the game G on tt^, and is a solitaire game whose value can be easily computed 
by means of a simple variant of Bellman-Ford algorithm. Therefore, the Ehrenfeucht and Mycielski’s 
results are already sufficient for determining a simple exponential time algorithm computing the node 
values and the two optimal positional strategies in an MPG: the algorithm consists in evaluating 
each possible strategy for one of the two players as a solitaire game for determining the optimal one. 
In the literature there are many local search algorithms that explore this space in a more efficient 
way [5,6,41,42] and some of them have been proven to be practically efficient in many settings by 
experiments [6,41]. Moreover, the global optimization problem of computing the best strategies for 
one player, according to a given metric, has been shown to have the property that every local optimum 
is also a global one for many complete metrics [5[. 

As another line of research, Zwick and Paterson [47] proposed pseudo-polynomial time algorithms 
for computing values of games, as well as positional optimal strategies. In particular, they considered 
the following four algorithmic problems: 

1. MPG-Decision(^, s): given a real number v and a start position s, decide whether Ug > v. 




14 



2. MPG-Threshold(T): given a real number T, determine for which nodes s G P it holds that 
Vs > T; 

3. MPG- Values: compute the optimal values Vg for all s G V. 

4. MPG-SynthesiS: assuming Vg > 0 (vg < 0) for every s G V, synthesize a positional winning 
strategy for Player 1 (Player 0); 

and they proved the following theorem: 

Theorem 5 ([47]). Let G = (V,E) be a mean payoff game. Assume all weights are integers and let 
W = maxg^E |w^(e)|- Then the following hold: 

1. MPG-Threshold(T) ean be solved in time 0(|Vp|i?| W) when T G Z, whereas it can be solved 
in time 0(|Vp|£^| W) when T G R; 

2. MPG- Values can be solved in time 0(|Vp|i?|lV); 

3. MPG-Synthesis can be solved in time 0(|V|'^|£^| log(|i?|/|V|) TV). 

Then, they observed that MPG-Decision is the basic decision problem for MPGs in the sense 
that several natural questions for MPGs, like evaluating the value Vg for every node s or constructing 
the optimal positional strategies, may all be Turing-reduced to it. They also pointed out that the 
existential results of Ehrenfeucht andMycielski [21] already implies that MPG-Decision G NPflco-NP 
and asked whether there might exist a strongly polynomial time decision procedure. Proving the 
existence of such algorithm is an open problem [7[. Finally, they showed how to reduce mean payoff 
games to other important families of games on graphs, like discounted payoff games and simple 
stochastic games. 

The complexity status of MPG-Decision has been since updated by proving that it lays in 
UP n CO-UP by Jurdzihski in [27[. 

In recent years, some other interesting results have been proven. Notably, in 2007 Lifshits, 
Pavlov [33] proposed a potential theory for MPGs and in 2011 Brim et al. [7] obtained faster al¬ 
gorithms exploiting results obtained in the the fields of Energy Games and energy progress measures, 
which are intimately related to the potentials studied in [33]. 

Their algorithmic results are summarized in the following theorem. 

Theorem 6 ([7]). For MPGs in which all weights are integers and for T G Z, the Value Iteration 
Algorithm [7] solves MPG-Threshold(T) and MPG-Synthesis in time 0(|P[ \E\ W), where W = 
maxegs |w^(e)[- 

We remark that both the algorithm of Paterson and Zwick [47] and the Value Iteration Algo¬ 
rithm [7] prescribe well defined procedures even if the weights on the arcs are real values. What is 
lost in running these algorithms on real weights is only the pseudo-polynomial upper bound on their 
running time. 

For our purposes, the family of pseudo-polynomial algorithms for MPGs is the best option. In¬ 
deed, in most of temporal workflow graphs all weights are expressed by integers of relatively small 
magnitude with respect to the intrinsic temporal granularity of the considered workflow. For exam¬ 
ple, in a workflow containing temporal distance constraints of days, the commonly adopted temporal 
granularity is the “minute” (m) and, therefore, all weights can be assumed to be less than 10^ as 
order of magnitude. In such circumstances. Brim’s algorithm offers the guarantee to terminate within 
short computation times. For these reasons we opted for integrating the procedures of Zwick and 
Paterson, as well as the faster procedures of Brim et al. [7], in order to efficiently solve instances of 
HyTN-Consistency and compute feasible schedules. 

Furthermore, as will be discovered in the experimental section, if these algorithms are suitably 
adapted—so as to allow them to terminate earlier as soon as certain evidences of inconsistency have 
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been collected—then their observed behavior outperforms by orders of magnitude what predicted 
by their theoretical pseudo-polynomial bounds even on input instances containing very large integer 
values. 

Based on these findings, we think that these pseudo-polynomial algorithms are to be considered 
(and probably adopted) even for solving HyTN instances where weights are floating point values whose 
magnitudes may differ in a significant way. In case the running time results to be unacceptable for a 
real application, one could then consider the possibility to round the weights to integer values. This 
rounding would clearly require special care: a very accurate approximation might lead to very high 
computation times while a gross approximation might not represent the original instance in a correct 
way. 


6 The Reductions 

This section presents the direct connection and the computational equivalence between MPG-Threshold 
and HyTN-Consistency. The equivalence is formally proven by offering one reduction in each di¬ 
rection. 

The reduction of HyTN-Consistency to MPG-Threshold allows to apply, in the context of 
HyTNs, any of the algorithms known for MPGs, included the exponential and subexponential ones. 

7 

Vice versa, in consideration of the fact that the MPG-Decision G P question is an open prob¬ 
lem [5,7,27,41,47], the reduction of MPG-Decision to HyTN-Consistency confirms that HyTN- 
Gonsistency offers an algorithmically more ambitious and mathematically steeper generalization of 
STN-Consistency (see also Remark 1). Moreover, the reduction gives a further evidence that, within 
STNs, a new algorithmic approach is necessary in order to manage temporal aspects of event like the 
synchronization one presented in the Introduction. 

Let us start considering the first reduction. 

Theorem 7. There exists a log-space^, linear-time, local-replacemeniP reduction from HyTN-Consistency 
to MPG-Threshold. 

Since this reduction plays a main role in the algorithmic solutions proposed in this paper, we firstly 
describe how it works and, secondly, we prove its correctness by means of two lemmas. Lemma 4 and 
Lemma 5. 

The reduction goes as follows. 

Let TL = {V, A) be a HyTN. We assume that every ri G V is the tail of some arc A G A. This 
assumption is not a restriction since, if "H contains a sink node v, i.e., a node v with no arc A G A 
having tail in it, then T-L is consistent if and only if so is Hy, the HyTN obtained from T-L by removing 
node V and every hyperarc having v as an head. Indeed, any feasible scheduling s : V i—>■ R for TL, 
once projected onto V \ {u}, gives a feasible scheduling for Tiv since every constraint involving v has 
been dropped and no constraint has been added; conversely, any feasible scheduling s for TLv can be 
easily extended to a feasible scheduling s for TL by exploiting the property of v being a sink node: it 
is sufficient to set s{v) = min{s(t/i) — wa{v) \ A G A, v G Ha}- 

Now, let us consider a mean payoff game Gy, = {Vo AVi,E) where: (1) Vo = V, Vi = A, nodes in 
Vo are colored by black while nodes in Vi are colored by white, and (2) for each A G A, the following 
weighted arcs are added to E: 

• an arc of weight 0 from the black node ta to the white node A, i.e., arc {tA, A, 0); 

• for each head node h G Ha, an arc of weight WA{h) from the white node A to the black node h, 
i.e., arc {A, h,WAih)). 

strong and basic-form of reduction introduced by Papadimitriou in |38]. 
restricted kind of Karp reduction introduced in |24]. 
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(a) Hyperarc A. 


(b) MPG representation of A. 


Figure 7: The conversion of a hyperarc into a white MPG node and its incident arcs. 


Algorithm 1: makeACorrespondingGcmie(?^) 
// a HyTN n = (V, A) 

1 

2 ^ A; 

3 £^ 0 ; 

4 foreach A € A do 

5 E <— E U {tA, A,0)-, 

6 foreach h € Ha do 

7 EE U {A,h,WA{h))-, 

Output: The MPG Gn = (Vb G Vi,E) 


Figure 8: The algorithm implementing the reduction from a HyTN to the corresponding MPG. 


In short, G-h = {Vq U Vi, A), with Vq = V, Vi = A, E = 0) | A e A} U {{A,h,WA{h)) \ 

A € A, h G Ha}- Fig. 7 depicts how a hyperarc is transformed into a MPG subnetwork while Fig. 8 
reports a pseudocode for the whole construction process, i.e.. Algorithm 1. 

G-u has \V\ + |A| nodes and 0{m) arcs and can be constructed in linear time. Moreover, G-h is a 
bipartite graph with bipartition (Vq, Vi) and it has been obtained from 77 by a simple local replacement 
rule: replace every hyperarc A G A by a claw subgraph as depicted in Fig. 7. For each single object, it 
is necessary only to manage a constant number of indexes, each of them having a polynomial size; thus 
the reduction is log-space. Fig. 9 depicts an MPG obtained applying the reduction to the motivating 
example HyTN depicted in Fig. 3; 

Now, let us introduce the formal proof of Theorem 7 by the following two lemmas. 

Lemma 4. If H is consistent then every node of G-h is a winning start position for Player 1 . 

Proof. Since H is consistent, there exists a feasible scheduling s : V —>■ R such that, for each hyperarc 
A G A, the reduced slack weight is non-negative > 0. Gonsider the following positional strategy 
TTi for Player 1: for each A G Vi, 

7ri(A) = arg min {s(h) — WA(h)}. 

IiGHa 

We claim that tti ensures Player 1 the win, wherever node the game starts from and however Player 0 
moves. In order to show this, we prove that the projection graph Gtti is conservative exhibiting a 
feasible potential p. Let p : Vo U Vi —>■ R be defined as follows: 


p{v) = 


s{v) 

s{t{v)) 


if « e Vo, 
if r: G Vi. 


Now, let a = {u, v, w) be any arc of Gt^ ■ 


( 2 ) 
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Figure 9: The MPG equivalent to the HyTN depicted in Fig. 3. A winning positional strategy tti for Player 1 
is highlighted by thick arcs. The dashed arcs are those not prescribed by strategy tti, i.e., they are removed 
when projecting the MPG on tti. 

Case 1: ii v G Vi, then w is a hyperarc of T-L with t{v) = u and w = 0; therefore, p{v) = s{t{v)) = 
s{u) = p{u) since u G Vq- Then wP{u,v) = w — p{v) + p{u) = 0 > 0 follows; 

Case 2: if S Vq, then u G Vi and w = Wu{v). Moreover, v = tti^u), which implies that v = 
argmin^gj:/^{s(/i) — w„(h)}. Therefore, recalling that > 0, i.e., s{t{u)) > min/ig//„{s(/i) —Wu{h)}: 

p{u) = > min {s(/i) — Wu{h)} = s(v) — Wu(v) = p(v) — w. 

hGHu 

Hence, w'^{u,v) = w — p(v) + p(u) > 0. 

In conclusion, G,ri is conservative. Therefore, the positional strategy tti certifies that any node of G 
is a winning start position for Player 1. □ 

Lemma 5. If every node of G-h is a winning start position for Player 1 then P is a consistent HyTN. 

Proof. If every node is a winning start position for Player 1, then there exists a positional strategy 
TTi which is everywhere winning for Player 1. Notice that must be conservative since Player 0 
can clearly win any play starting from a node located on a negative cycle. Let p : Vq U Hi —> R be 
a feasible potential for G,ri • We claim that the restriction of p onto Vq is a feasible scheduling for 
H. Indeed, for any hyperarc A of H, {tA,A,0) is an arc of Gtti, whence p{A) < p{tA)- Moreover, 
(H,7ri(H),'u;^(7ri(H))) is also an arc of G,ri, whencep(7ri(A)) < p{A)+WAi'n'i{A)). Since7ri(H) G Ha, 
then the following holds: 

pitA) > p{A) > p{Tri{A)) - WAiTTi{A)) 

> min {s(/i) — 
hGHA 

Hence, the restriction of p onto Vq is a feasible scheduling for TL. Thus, TL is consistent. □ 

In Fig. 10 the values under the nodes represent a feasible potential for the projection of the MPG 
depicted in Fig. 9. By Lemma 5, the restriction of such a feasible potential on the black nodes is also 
a feasible scheduling for the corresponding HyTN depicted in Fig. 3. Now, we have all the necessary 
results to prove the following theorem. 
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Figure 10: The integer labels under the nodes are a feasible potential for the projection on tti of the MPG 
depicted in Fig. 9. The restriction of this potential on the black nodes (those in Vb) is a feasible scheduling 
for the HyTN depicted in Fig. 3 as explained in the proof of Lemma 5. 


Algorithm 2: isConsistent(?^) 

//a HyTN H = (V,A) of unknown consistency state 

1 Gn <—makeACorrespondingGame (Lf); // See Algorithm 1 

2 (VFojlFi) -f- solveMPG-THRESHOLD(G'H, 0); // Brim’s algorithm, see Theorem 6 

3 if (Wo = 0) then Output: YES; 

4 else Output: NO; 


Figure 11: Pseudocode of the algorithm for deciding HyTN-Consistency. 


Theorem 8. LetTi. = {V,A) he an integral-weighted HyTN, m = o,ndW = |tyA(^)|} 

the maximal weight value present in TL. The following propositions hold: 

1. There exists an 0((|y | + |A|)mLF) pseudo-polynomial time algorithm deciding HyTN-Consistency 
for H; 

2. There exists an 0((|y| + |A|)mLF) pseudo-polynomial time algorithm such that, given on input 
any consistent HyTN TL, it returns as output a feasible scheduling s : Vu -A- Ti of TL; 

3. There exists an 0((|y| + |A|)mLF) pseudo-polynomial time algorithm such that, given on input 
any not-consistent HyTN TL, it returns as output a negative cycle (S,C) ofTL. 

Proof. 1. The decision algorithm is sketched in Fig. 11. It takes in input a HyTN TL = {V,A) and, 
in line 1, constructs the corresponding MPG Gn as described in Theorem 7. This first step 
takes 0(m) time and yields a graph with \V\ + |A| nodes and 0(m) arcs. Then, in line 2, the 
instance of MPG-Threshold with T = 0 on graph G-h is solved in 0((|y | + |A|)to 1T) time by 
the Value Iteration Algorithm (see Theorem 6). The output consists in a partition of Gn nodes 
into two sets: Wi = {v € V U A \ Vv > Q} and Wq = {u G V U A | < 0}. If Wq is empty, then 

TL is consistent by Lemma 5, otherwise it is not consistent by Lemma 4. 
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Algorithm 3: computeAFeasibleSchedule('H) 

//a consistent HyTN H — {V, A) 

1 GmakeACorrespondingGaine(?^); // See Algorithm 1 

2 TTi MPG-Synthesis (G); // Compute a positional winning strategy for Player 1; see Theorem 6 

3 Gtti t— compute the subgraph of G induced by tti; 

// Recall G^i = {Vo U Vi, E), where Vo = V and Vi = A. 

4 s a new node; // s ^ Vo UVi 

5 Add s to Vi and add an arc (s, v, 0) for each v £ Vo', 

6 p ■£- Bellman-Ford(G 7 ri, s); // compute a potential function p 
Output: the restriction of p onto V 


Figure 12: Pseudocode of the algorithm for computing a feasible schedule. 


2. In case Wq is empty, a feasible scheduling is obtained as shown in Algorithm 3. First, in line 2, the 
algorithm computes a positional winning strategy tti for Player 1. This takes 0((|P| + |A|)toIT) 
time by Theorem 6. Next, in line 3, it builds the graph which is conservative since tti is a 
positional winning strategy for Player 1. Then, in lines 4-5, it adds a new node s to Vi and a 
new arc = (s, v, 0) for each node v £Vo\n . Let = {Vq U (Vi U {s}), E') the graph thus 
obtained. Observe that every node of G'^_^ is reachable from s. Indeed, every node A £ Vi = A 
can be reached by traversing two arcs: from s to Ia along the arc = (s, Ia, 0), which belongs 
to G(rj as Ia £ Vo, then from tA to A along the arc {tA, A, 0), which belongs to G^^ (and hence 
to G(rJ since G Vq. 

Since the added node s is a source, then G'^^ is conservative too. Therefore, in G'^^, the set of 
distances from node s, computed calling the Bellman-Ford algorithm in line 6, forms a feasible 
potential p : Vq U Vi U {s} ^ Z and the restriction of p onto Vq = V^ is a feasible scheduling for 

n. 

3. In case Wo is not empty, a negative cycle is determined by Algorithm 4. Let G[VVo] be the 
subgraph of G induced by VTo) be., the graph obtained from G by removing all nodes not in 
VTo and all the arcs incident into them. Notice that every node v £ Wo is a winning start 
position for Player 0 in game G[VVo] because u is a winning start position for Player 0 in game 
G, and no winning strategy for Player 0 in G can prescribe a move from a node in VVo to a 
node in VFi; therefore, that same winning strategy remains valid on G[Wb]. This implies that, 
for every u £ Wo, there exists at least one arc {u,v) with v £ Wq- In particular, since {Vo,Vi) 
is a bipartition of G, then Wo = Wq fl V/q ^ 0. In line 3, a positional winning strategy ttq for 
Player 0 on G[VFo] is determined. By Theorem 6, this computation takes time 0((|V^|-|-|A|)mVT). 
Consider the set of hyperarcs C = {7i'o(f)}^gn7Q; the pair (Wo,C) returned by the algorithm is 
a negative cycle. Indeed, for any v £ VTo) '^o(v) € Vi is a hyperarc of Thus the head set 

^ Vq. Also, C VTo, since u is a winning start position for Player 0 and ttq is a 

winning strategy for Player 0. Combining, C VFq determining that {Wo,C) is a negative 

cycle. 

□ 

Remark 2. In Theorem 8 Item 2), a set of feasible potentials may be obtained without executing 
the Bellman-Ford algorithm. Actually, if the partition {Wo,Wi) is computed by the Value Iteration 
Algorithm [7], then a feasible scheduling for "H can be directly derived from the progress measure 
computed within the algorithm. In more detail, let G = {Vq U Vi,E) be an MPG weighted by 
w : E ^ Z. An energy progress measure is a function / : Vq U V4 —t N U {-foo} such that: if 
V £ Vo, then for every {v,v',w) £ E it holds f{v) > f{v') — w, otherwise, v £ Vi and there exists 
{v, v', w) £ E such that f{v) > f{v') — w. An energy progress measure / : Vq U V4 ^ N U {+cxd} such 


20 





Algorithm 4: computeANegativeCycle('H, M/q) 

//a HyTN 'H = {V,A) = (Vb U Vi,A) which is not consistent 
// the non-empty set Wo — {v G V \ Vy < 0} 

1 G ■<—inakeACorrespondingGaine(?^); // See Aigorithm 1 

2 G[Wo] t— compute the subraph of G induced by Wo; 

3 no G- MPG-Synthesis (G[VKo]); // Compute a positional winning strategy for Player 0; see Theorem 6 

4 W 0 G- Wo n Vb; 

5 C <r- {rro(u)}^g-jy^; 

Output: (WoyC) 


Figure 13: Pseudocode of the algorithm for computing a negative cycle. 


Algorithm 5: computeAFeasibleSchedule-Remark2(iA) 

//a consistent HyTN H — (P, A) = (Po U Pi, A) 

// ref. Remark 2 and Theorem 6 [7] 

1 G makeACorrespondingGaine(?f); // ref. Algorithm 1 

2 / -f- Value-Iteration(G); // compute an energy progress measure for G as in Theorem 6 
Output: / 


Figure 14: Pseudocode of the algorithm of Remark 2 for computing a feasible schedule. 


that 0 < f(v) < -boo for every t; C Pn U Pi is provided by the resolution algorithm of Theorem 6 in 
time 0((|P| -b |A|)to1P). 

The progress measure / is already a feasible scheduling for H: in fact, for every hyperarc A C A, 
it holds {tA,A,0) € E and {A,v,wa{v)) G E, for every v € El a; combining these two last facts, it 
follows that: 

f{tA) > f{A) > min {/(u) - u;A(t^)}, 
vGHa 

i.e., / is a scheduling satisfying all constrains A € A. This allow us to employ the algorithm depicted 
in Fig. 14 instead of the one depicted in Fig. 12 in the case that IPi = P. 

The computational equivalence between MPG-Decision problem and HyTN-Consistency can 
be now determined by showing that also MPG-Decision can be reduced to HyTN-Consistency. 

Theorem 9. There exists a\og-space, linear-time, local-replacement reduction from MPG-Decision 
to HyTN-Consistency. 

Proof. Let G = (Pq U Pi,i?) be an MPG. For each node w € Vb U Pij let Na{u) denote the outgoing 
neighborhood of u in G, i.e., Nc{u) = {u G Vb U Pi | (u, u) G E}. 

A corresponding HyTN H = (P, A), where P = Pq U Pi, is constructed from G as follows. For 
every u G Pi, a hyperarc A„ G A is added to TL, where: 

A„ = {u,Ng{u),waJ, 

with weight wAu{v) — w{u,v) for every v G Ng{u). Moreover, for every u G Vb and every v G Ng{u), 
a hyperarc Auv G A is added to TL, where: 

Auv = iu,v,w{u,v)). 

This construction requires a log-space and linear-time computation. 

Now, we firstly prove that if TL is consistent then every node of G is a winning start position for 
Player 1. 
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Indeed, let s : F —s- R be a feasible scheduling for T-L. Thus, w\ > 0 for every hyperarc A ^ A. 
Notice that, by construction, for each u &Vi there exists a unique hyperarc € A with tail = w; 
moreover, it holds that Ha^ = Ng{u). Hence, for each u € Vi, we can define a positional strategy tti 
for Player 1 as follows: 

7ri(u) = arg min {s{h) — WA^ih)}. 
h&HAu 


Now, consider the potential function p : P —s- R defined as: p{u) = s(u) for every u € V. We argue 
that p is a feasible potential for ■ 

In fact, let a = (u, v,w) G E be any arc of G: 


Case 1: Assume that u € Vq. Then, by construction, a = A^v- Hence, from p 
of s, we have: 


p(u) = s{u) > minheHAuvi^i^) ~ 
= s{v) - WA^^iv) 

= p(v) — w 


A 


s and the feasibility 


Hence, w'p(u, v) = w — p(v) + p(u) > 0; 


Case 2: Assume that u G Pi. Then, by construction, A = (u, Ng(u),wa), where wa{z) = w{u,z) for 
every z G Ng{u)] moreover, notice that v = 'Ki{u) G Ng{u) = Ha- Hence, from p = s, the 
feasibility of s, and the definition of tti, we have: 


p{u) 


s{u) > - WA^{h)} 

= s(7ri (m)) - WA ^ (tti (u)) 

= s{v) -WA^iv) 

= p{v) — w 


Hence, w'p(u, v) = w — p(v) — p(u) > 0. 


Thus, Ctti is conservative. This implies that every node of G is a winning start for Player 1. 

Secondly, we prove that if every node of G is a winning start position for Player 1, then H is 
consistent. 

Let TTi be a positional winning strategy for Player 1. It follows that G,ri is conservative and, 
therefore, it admits a feasible potential p : P —s- R. Now, consider the scheduling function s : P —>■ R 
for H defined as: s{u) = p(u) for every u G V. We argue that s is a feasible scheduling of H.. 

In fact, let A= {Ia, Ha, wa) G A be any hyperarc of H: 

Case 1: assume Ia G Vo- Then, by construction, A = (u, v, w) for some v G Ng{u),w G R and u = tA- 
Hence, from s = p and the feasibility of p, we have: 


s{tA)=piu) > p{v)-W 

= s{v)-WAiv) 

= min/,g_f/^{s(h) - wa(/i)} 

Hence, s satisfies A, i.e., w%>0 -, 

Case 2: assume Ia G Pi. Then, by construction, A = {u, Ng{u),wa) for u = tA and wa{v) = 
w{u,v) G R for every v G Ng{u)', moreover, if u = 7ri(t6), then v G Ng{u) = Ha- Hence, from 
s = p and the feasibility of p, we have: 

s{tA)=piu) > p(v)-w 

= s{v)-WAiv) 

> min/,g_f/^{s(/i) - WAih)} 

Hence, s satisfies A, i.e., u;^>0. 

This proves that s satisfies every hyperarc A G A- Then s is a feasible scheduling of H, which is thus 
consistent. □ 
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7 Computational Experiments 

This section describes our empirical evaluation of the proposed consistency checking algorithms to 
evaluate the performances and the general applicability of the proposed HyTN model. Both Algo¬ 
rithm 3 and Algorithm 4 consist of one single call to Algorithm 2, plus some extra computation 
of lower asymptotic complexity. Since the cost of these further computations was confirmed to be 
practically negligible in some preliminary experiments, we report on the results of our experimental 
investigations only for Algorithm 2. 

All algorithms and procedures employed in this empirical evaluation have been implemented in 
C/C++ and executed on a Linux machine having the following characteristics: 

• 2 CPU AMD Opteron 4334; 

• 64GB RAM; 

• Ubuntu server 14.04.1 Operating System. 

The source code and all HyTNs used in the experiments are freely available [14]. 

The main goal of this empirical evaluation was to determine the average computation time of 
Algorithm 2, with respect to randomly-generated HyTNs following different criteria, in order to give 
an idea of the practical behavior of the algorithm. According to Theorem 8, the worst-case time 
complexity of Algorithm 2 is 0{{n + m')mW), where n = \V\, m' = |A|, m = 

W = max^g_ 4 {max/ig^ |w^(ft,)|}. Hence, we implemented different experiments with respect to the 
parameters and W. Here we propose a summary of the obtained results presenting a brief 

report about four tests. Test 1, Test 2, Test 3 and Test 4. 

In Test 1 the average computation time was determined for different HyTN orders n to empha¬ 
size the practical computation time dependency on n. In Test 2 the average computation time was 
determined for different HyTN maximal edge-weights W to understand how much the practical com¬ 
putation time is dependent on W. In Test 3 we investigated how some execution times affect the 
value of the standard deviation, with the goal to determine how many instances require a significant 
greater computation time with respect to the average time of a data set. Finally, in Test 4 the average 
computation time was determined with respect to different values of the number of possible strategies 
of Player 1 0/16+ order to give an idea about the possible practical relation between execution 

time and number of possible strategies. 

The generation of random HyTN instances was carried out exploiting two generators. The first 
generator was the random workflow schema generator provided by ATAPIS toolset [31]: it produces 
random workflow graphs according to different input parameters that allow to control the minimal 
and maximal number of activities, probability for having parallel branches, the minimal and maximal 
probability of inter-task temporal constraints, etc. on the generated graphs. We verified that this tool 
allows the determination of graphs that are not only a closer approximation to real-world examples, 
but also more difficult to check than those generated at random without particular criteria. 

We generated benchmarks as follows: 

1. First, temporal workflow graphs were generated by fixing the probability for parallel branches 
to 10% and maximal value for each activity duration or delay between activities to a value W, 
where W was chosen accordingly to the test type; 

2. Then, each workflow graph was translated into an equivalent HyTN Ti. by the simple transfor¬ 
mation algorithm exemplified in Section 2. 

It is worth noting that different random workflow graphs all having the same number of activities may 
translate into HyTNs having different orders n because the original workflow graphs may have different 
number of connector nodes. Considering the transformation algorithm exemplified in Section 2, it is 
easy to verify that a workflow with N activities can translate into a CSTN having between 2N + 2 
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Table 1: Orders of the smallest and biggest HyTN determined for each set of random generated workflows 
having N activities. 


N 

Order of smallest HyTN 

Order of biggest HyTN 

10 

26 

50 

20 

48 

94 

30 

78 

138 

40 

104 

196 

50 

136 

236 

60 

164 

268 

70 

196 

306 

80 

222 

350 

90 

262 

394 

100 

292 

410 


Table 2: Comparison between different kinds of queue implementation in the Value-Iteration procedure. All 
values are in seconds. 



FIFO Queue 

LIFO Queue 

LIFO Queue 
+ Stopping-Griterion 

Max-Priority Queue 


90.55 

11.77 

6.98 

184.69 

a 

487.69 

64.10 

34.61 

653.26 


nodes (when the workflow is a simple sequence) and 5N -k 2 nodes (when the workflow is a sequence 
of groups of two parallel activities). 

ATAPIS toolset has been designed to generate graphs with strongly connected components (An¬ 
dreas Lanz, personal communication, October 6, 2015). In particular, it has been optimized for small 
graphs with up to 50 activities. This design choice was motivated by the widely accepted seven process 
modeling guidelines [34] which suggests to always “decompose a (workflow) model with more than 50 
elements (activities)”. Therefore, we used the tool for generating random workflow graphs with 100 
activities at maximum and, consequently, obtaining HyTNs having 502 nodes at most. 

In Table 1 we report the orders of the smallest and the largest HyTN determined from each set of 
random generated workflow graphs having N activities for N G {10, 20,..., 100}. 

In order to study the scalability of the algorithm with respect to the number of nodes, we had 
to rely on a second generator of random HyTN graphs. Our choice has been to use the randomgame 
procedure of pgsolver suite [39], that can produce parity games instances for any given number of 
nodes. In particular, we exploited randorngame in the following way: 

1. First, randorngame was used to generate random directed graphs with out-degree fixed to 3; 

2. Then, the resulting graphs were translated into MPGs by weighting each arc with an integer 
randomly chosen in the interval [—W, W], where W was chosen accordingly to the test type; 

3. Finally, each MPG G was translated into a HyTN TLg by the reduction algorithm of Theorem 9. 
During the translation from MPG to HyTN, only 10% of the hyperarcs were maintained having 
multiple heads, while 90% of hyperarcs were transformed into standard arcs. This 10%-rule 
stems from the fact that we are considering workflow based applications where the percentage 
of (multi-headed) hyperarcs is less than 10% compared to standard arcs in general. 

With such settings, the resulting HyTNs are characterized by m,m' G 0(n). 

Before presenting the summary of results, it is worthwhile to present some implementation choices 
about Algorithm 2 that we had to adopt. The core of the algorithm consists of calls to algo¬ 
rithms makeACorrespondingGame(7^), that transforms the given HyTN % into a MPG Gn, and 
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solveMPG-Threshold(G-H, 0) (Value Iteration algorithm), that determines for which game nodes s 
it holds that Vs > 0. The makeACorrespondingGame () implementation didn’t require significant 
choices thanks to the simple structure of the algorithm. On the contrary, in the implementation of 
solveMPG-ThresholdO we introduced some further ideas in order to speed-up the algorithm and 
avoid unnecessary computations. In particular, it is not necessary for solveMPG-ThresholdO to con¬ 
tinue to determine other potential value Vs' as soon as it determines a value Vs < 0: at this point we 
can already conclude that the network is not consistent and, with a lower computational cost, we can 
yield a generalized negative circuit assessing this fact (Lemmas 4 and 5). 

Moreover, we verified that there is an important data structure in the original Value Itera¬ 
tion algorithm, a queue, that is not further specified by the authors and that different implemen¬ 
tations of it affect the performance of the algorithm. Therefore, we decided to verify whether 
solveMPG-ThresholdO performance could be appreciably improved adding a suitable stopping cri¬ 
terion and a proper queue implementation. Table 2 reports the obtained results, mean execu¬ 
tion time fx and its standard deviations cr, determined running the following different versions of 
solveMPG-ThresholdO on the same data set of 10^ not consistent HyTNs^ having \V\ = 10® and 
W = 10®: 

1. FIFO Queue: the original queue is implemented as a FIFO queue; 

2. LIFO Queue: the original queue is implemented as a a LIFO queue (stack); 

3. LIFO Queue-|-Stopping Criterion: the queue is implemented as stack and the computation 
is halted either when all potential values are stable or when any of them is negative; 

4. Max-Priority Queue: the original queue is implemented as a Fibonacci’s heap. 

The results show that, in general, solveMPG-ThresholdO performance becomes better if the 
original queue is implemented as a stack and, in particular, a further improvement can be obtained if 
the stopping criterion is also considered. Nevertheless, such improvements can only partially reduce 
the statistics variability of the running time, as it is shown in the following experimental results. 

As mentioned above, the goal of Test 1 was to determine the average computation time of Algo¬ 
rithm 2 implementation for different values of n to study the practical computation time dependency 
on such parameter. 

The instances in Test 1 come from the remdorngcune generator, except those for the first row of the 
table in Fig. 15a which have been built by the ATAPIS workflow random generator. In particular, 
for each n G {1 • 10®, 2 ■ 10®,..., 10 • 10®}, 1600 HyTN instances with maximum weight W := 1000 
and unknown consistency state were generated by randomgame, whereas 1600 HyTNs of unknown 
consistency state and order n around 400 were generated by ATAPIS. The results of the test are 
summarized in Fig. 15, where each execution mean time is depicted as a point with a vertical bar 
representing its confidence interval determined according to its standard deviation. 

The depicted function interpolating the mean values shows that the practical performance of the 
algorithms is definitely better than the theoretical worst-case bound of 0((n + m')mW); in our case 
this last is 0{'n?) since in the generated data sets W is constant and m,m' € 0(n). Fig. 15c depicts 
the interpolating function of experimental execution times and, in red, the function n^/10^® as a 
reasonable surrogate for the worst-case execution time. The comparison shows that the algorithm 
performs very well in real case executions. 

However, since the standard deviation observed in the experiment is not negligible, we further 
investigated the behavior of the algorithm and we discovered that there is a correlation between the 
execution time of the algorithm and the consistency state of the input HyTN. Therefore, /i and a 
were recalculated considering two kind of HyTN sets: one having all consistent HyTNs, and the other 
having all not consistent HyTNs. 

"^We considered not consistent HyTNs because they practically required more time to be solved. 
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n 

/i (sec) 

a 

<4-10^ 

0.13 

0.42 

1 • 10® 

0.55 

5.41 

2-10® 

0.99 

4.71 

3-10® 

1.67 

13.55 

4-10® 

1.95 

12.59 

5-10® 

2.58 

16.10 

6-10® 

2.58 

9.43 

7-10® 

3.48 

22.43 

8-10® 

4.58 

17.85 

9-10® 

4.72 

36.19 

10 • 10® 

4.83 

30.62 



(a) Test 1 results. (b) Interpolation of average execution times 

of Test 1. 



(c) Comparison between theoretical computation times and experimental ones. 


Figure 15: Results of Test 1: average execution times (^) and relative standard deviations (a) over a range of 
different HyTN orders n. Times are in seconds. Each data set comprised of 1600 HyTN instances of unknown 
consistency state. 


Fig. 16 depicts average execution times obtained in Test 1 calculated considering samples of either 
all consistent or all not consistent HyTNs obtained from workflow graphs. Fig. 17 offers the same 
view but for HyTNs obtained from MPG graphs. In general, the mean execution times for consistent 
HyTNs are smaller than the corresponding ones for not consistent HyTNs; furthermore, they also 
exhibit a negligible standard deviation. However, for samples of consistent HyTNs obtained from 
workflows, the standard deviation is not negligible even for samples with size N = 20. Part of the 
reasons for this behavior is given by the structure of the data sets: in each data set HyTNs can differ 
a lot with respect to their order and, therefore, they may require very different execution times. For 
example, the data set relating to workflow graphs with 20 activities contains HyTNs with order in 
range [48,94]. Since the number of activities is usually considered as main parameter in workflow 
community, we wanted to maintain such structure of data set and experiment results to emphasize 
the dependency of execution time with respect to such number. 

On the other side, for consistent HyTNs determined from MPGs, the observed standard deviation 
a is always less that the 10% of the average execution time /r with 99% level of confidence, while for 
not consistent HyTNs it has not been possible to determine any confidence level because the observed 
standard deviation a resulted to be always high due to some hard instances. 

Even though procedure solveMPG-ThresholdO could require up to Q{W) updates according to 
the theoretical worst-case bound, our experiments suggest that, in practice, some kind of dependency 
of the running time on W is appreciable only for a few MPG games, all associated to not consistent 
HyTN instances. 

The goal of Test 2 was to determine the average computation time of Algorithm 2 for different 
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N 

/i (sec) 

a 

10 

6.42 • 10“'' 

1.22 • 10“'’ 

20 

1.05 • lO-'^ 

4.85 • 10“® 

30 

1.50- lO-'^ 

5.7-10“® 

40 

2.43 • lO-'^ 

1.04 • 10“® 

50 

3.20 • lO-'^ 

1.78 • 10“® 

60 

3.77- lO-'^ 

1.38 • 10“® 

70 

4.77 • lO-'^ 

1.28 • 10“® 

80 

5.73 • 10“® 

1.80 • 10“® 

90 

6.82 • lO-'^ 

2.79 • 10“® 

100 

8.89 • lO-'^ 

4.10 • 10“® 



(a) Average execution times for consistent HyTfSi^ Interpolation of average execution times of Ta- 
obtained from workflow graphs with N activitiMe 16a. 


N 

fi (sec) 

a 

10 

4.45 • 10“® 

1.38 • 10“® 

20 

1.50 • 10“® 

5.10 • 10“® 

30 

4.04 • 10“® 

1.48 • 10“^^ 

40 

1.10 • 10“^^ 

3.62 • 10“^^ 

50 

1.64 • 10“^^ 

8.42 • 10“^ 

60 

4.36 • 10“^^ 

1.20 • 10“® 

70 

8.08 • 10“^^ 

2.71 • 10“® 

80 

1.31 • 10“® 

4.20 • 10“® 

90 

1.59 • 10“® 

5.22 • 10“® 

100 

2.59 • 10“® 

8.46 • 10“® 



(c) Average execution times for not consis^)i1Interpolation of average execution times of Ta- 
HyTNs obtained from random workflow gri^dflc. 
with N activities. 

Figure 16: Average execution times obtained in Test 1 calculated considering samples of either all consistent 
or all not consistent HyTNs obtained from workflow graphs. 


values of W, in order to understand how much the practical computation time is dependent on W. 
Therefore, we considered three possible edge weight ranges, [10^,10^], [10^,10®], and [10®, 10®], and 
for each of them two data sets have been built using the rcmdomgcune generator, one comprising only 
consistent HyTNs, and the other only not consistent ones. Each data set comprised of 800 HyTNs 
instances having \V\ = 10® nodes, m,m' € Q{n) and edge weights in the corresponding weight range. 
Fig. 18 depicts the results on these six data sets. Applying the worst-case analysis for these data sets, 
it results that the time complexity should be 0{W) since n, m and m' are constants. On the contrary, 
the determined interpolation functions representing the experimental execution times do not show any 
clear dependence on W. This result suggests that, in practice, uniform random weighted instances 
are decided very quickly with respect to the magnitude of their weights and that the algorithm does 
not seem to exhibit the worst-case pseudo-polynomial behavior predicted in the asymptotic analysis. 
Moreover, the average execution times for each data set comprising only consistent HyTNs are less 
than those for the corresponding data of only not consistent HyTNs. Only for consistent HyTNs data 
sets the standard deviation was below 7% than the average execution time with a confidence of 99%. 

In order to better understand how some execution times affect the value of the standard deviation, 
we conducted a third experiment. Test 3, with the goal to visualize the distribution of the instances 
with computation times significantly above the average. Procedure solveMPG-ThresholdO has been 
executed on 10® randomly generated not consistent HyTNs, each having order n = 10® and W w 10®. 
The determined running times are depicted in Fig. 19a: most of the instances are decided very quickly, 
i.e., in a time between 0 and 10 seconds, while a smaller portion of the HyTNs required a time between 
10 and 500 seconds. In more details, in repeated tests we verified that, approximately, 1% of the HyTN 
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(a) Average execution times for c0lilsibltelpolation of average execution times of 
HyTNs obtained from MPGs. Table 17a. 
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(c) Average execution times for n^4) clBfeigB@A»tion of average execution times of 
HyTNs obtained from MPGs. Table 17c. 


Figure 17: Average execution times obtained in Test 1 calculated for samples of either all consistent or all not 
consistent HyTNs obtained from MPG graphs. 


instances required a time between 50 and 100 seconds to be decided, 0.4% required a time between 
100 and 500 seconds, and, finally, only 0.1% required more than 500 seconds. These results are shown 
in Fig. 19b. 

Such behavior has been confirmed in other tests with different graph orders and maximum edge 
weight values. In several experiments we conducted, we observed that the maximum execution time 
of the algorithm keeps growing as we enlarge the size of the dataset. This explains why the standard 
deviation can’t be reduced. If we could characterize such hard instances in general, we would be 
making a major progress in understanding the computational complexity of MPGs. We didn’t find 
any pattern or property that characterizes the found hard instance. Here we can only show a simple 
family of HyTNs instances in which the execution time grows linearly with W. The family is given 
by just one single HyTN graph where only W changes, as depicted in Fig. 20a. The corresponding 
MPG is shown in Fig. 20b and provides a clear example where Brim’s Value Iteration algorithm [7] 
performs poorly. It is worth noting that in the context of MPGs this example can be reduced down 
to 6 nodes. 

Finally, in order to show how much the running time is dependent on the number of different 
positional strategies of one player, in Test 4 the average computation time has been calculated with 
respect to different values of the product of the heads of hyperarcs (i-e-, UagA\Ha\) in a HyTN. In 
particular, for each H G {10^^, 10^°,..., 10®^}, 2500 HyTNs instances {V,A) each having \V\ = 50 
nodes, ~ H, and W = 10® have been generated by means of randorngame generator. 

The results of the evaluation are depicted in Fig. 21, where H values are drawn in logarithmic scale. 
Analyzing the diagram in the figure it is possible to say that, experimentally, the average execution 
time increases only logarithmically with respect to the number of different positional strategies of one 
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Figure 18: Average execution times in Test 1 calculated considering samples of either all consistent or all not 
consistent HyTNs. 


player. This results is quite interesting because, considering the HyTN in Fig. 20a, it is evident that 
the time for checking a HyTN is more dependent on the edge weight magnitude than on the number 
of different positional strategies of one player. 

8 Related Work 

In the literature there are some extension proposals of the STN model to augment the capability to 
represent temporal constraints. 

In the STN seminal paper [18], Dechter et al. firstly proposed to consider the Temporal Con¬ 
straint Satisfaction Problem (TCSP). A binary constraint in a TCSP is represented using a set 
of intervals rather than a single interval as in an STN. In particular, a binary constraint Cij = 
{[oi, 6 i], [a2,b2], ■ ■ ■, [a;, 6 ;]} between time points Xi and Xj represents the disjunction oi < Xj — Xi < 
6 i V 02 < Xj — Xi < 62 V o; < Xj — Xi < bi. The problem of verifying consistency of a TCSP is 
NP-complete as the same authors showed in the paper; hence, they finally propose to consider STNs 
as a tractable simplified model. 

A similar kind of generalization considering disjunction of temporal distance constraints was 
proposed by Stergiou and Koubarakis [44] defining the Disjunctive Temporal Problem (DTP). A 
DTP consists of a set of variables X = {xi,X 2 ,... ,Xn} having continuous domains and represent¬ 
ing time points and a set of disjunctive difference constraints between the time points in the form: 
oi < xq - Xj^ < 61 V 02 < Xi^ — Xj^ < 62 V ... V Ofc < Xi^ — Xj^^ < bk\ where Xi^,Xj^,... ,Xi^,Xj^ 
are time points from X and ai,bi,... ,ak,bk are real numbers. A DTP is consistent if there exists an 
instantiation of variables X to real numbers satisfying all the constraints. Since DTPs are a general¬ 
ization of TCSPs, also for DTPs the consistency check problem is NP-complete. In [44] the authors 
presented some of the theoretical results characterizing the possible backtracking algorithms that solve 
the consistency problem in terms of search nodes visited and consistency checks performed. 

In 2005, Kumar proposed to consider a restricted class of DTP in order to maintain some of the 
expressive power of DTPs but, at the same time, allowing an efficient consistency check. In particular, 
in [40], RDTPs (restricted DTPs) is defined as a disjunctive temporal problem where a constraint is 
one of the following types: (Type 1) (/ < Xi—Xj < u), (Type 2) (/i < Xi < ui)V{l 2 < Xi < U 2 )---{lj < 
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(b) Instance classification with respect to 
solveMPG-ThresholdO execution time. 


Figure 19: solveMPG-ThresholdO execution times obtained in Test 3 determined considering samples of all 
not consistent HyTN instances. 


Xi < Uj), (Types) {h < Xi < ui) V (I2 < xj < U2), where Xi and xj represent a timepoint variable, 
and li,Ui real values. An RDTP instance can be solved in strongly polynomial-time deterministic 
algorithm transforming it into a binary Constraint Satisfiability Problem (CSP) over meta variables 
representing constraints of Type 2 or Type 3 and, then, showing that such binary constraints are also 
connected row-convex (CRC) constraints, and, then, exploiting the properties of CRC constraints. An 
instantiation of a consistency check algorithm for RDTPs that further exploits the structure of CRC 
constraints has a running time complexity of 0 ((|TP 2 | + l^-P 3 |)^c^max + + \TPz\Y{NM, 

where TP 2 is the set of Type 2 constraints, TP^ is the set of Type 3 ones, dmax is the maximum number 
of disjuncts in any constraint, and N/M is the number of the nodes/arcs of the instance, respectively. 
In the same paper, Kumar presented also a simpler and faster, but randomized, algorithm for the 
same class RDTP. 

An attempt to model some aspects of STNs similar to those addressed by HyTNs was lead in [2], 
where fun-in and fun-out subgraphs much resembling our multi-tail and multi-head hyperarcs were 
considered. However, since the problem l-in-3-SAT is NP-complete even when all the literals compris¬ 
ing the clauses are positive, it readily follows that their models lead to NP-complete problems even 
when fun-out subgraphs (or fun-in subgraphs) are banned. As such, the opportunity for tractability 
spotlighted in this paper is missed in those models. 

Another approach to extend STN is represented by the proposal of Khatib et al. [28,29]. They 
introduced the characterization of hard and soft constraints. STNs are able to model just hard 
temporal constraints, i.e., they can represent instances where all constraints have to be satisfied, and 
that the solutions of a constraint are all equally satisfying. However, such assumption can be too much 
restrictive in some real-life scenarios. For example, it may be that some solutions are preferred with 
respect to others and, hence, the main problem is to find a way to satisfy them optimally, according 
to the preferences specified. To address these kind of problems, in [28] the authors introduced a 
framework in which each temporal constraint is associated with a preference function specifying the 
preference for each distance or duration; a soft simple temporal constraint is a 4-tuple ((A, Y), I, A, f) 
consisting of (1) an ordered pair of variables (A, Y) over the integers, called the scope of the constraint; 
(2) an interval I = [a,b], where a and b are integers such that a < b; (3) a set of preferences A; (4) 
a preference function /, where / : [a, 6] 1 —> A is a mapping of the elements belonging to interval / 
into preference values, taken from set A. An assignment fa, and Vy to the variables A and Y is said 
to satisfy the constraint {{X,Y),I, A, f) if and only if a < — v^- In such a case, the preference 

associated to the assignment by the constraint is f{vy — Vx)- Using soft simple temporal constraint, 
a new model of temporal constraint network has been introduce: the Simple Temporal Problem with 
Preferences (STPP). In general, each solution of a STPP has a global preference value, obtained by 
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(a) A HyTN Hw for which Algorithm 2 takes <d{W) time. 



(b) The corresponding MPG computed by Algorithm 1. 

Figure 20: A HyTN which requires 0{W) computation time by Algorithm 2. 
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(a) Average execution times obtained for differ¬ 
ent values of the product of the heads of hyper¬ 
arcs n. 



(b) Interpolation of average execution times 
in Table 21a. 


Figure 21: Average Execution Times obtained in Test 4. 


combining in a suitable way the preference levels at which the solution satisfies the constraints. The 
optimal solutions of an STPP are those solutions which are not dominated by any other solution 
in terms of global preference. It was shown in [28] that, in general, STPPs belongs to the class of 
NP-hard problems. When the preference functions are semi-convex and some other side conditions 
are observed, then the problem to find an optimal solutions of an STPP is tractable [29]. 

Finally, another kind of possible extension is represented by the use of ^ operator instead of < 
in the binary constrains of STNs. Koubarakis [30] showed that if in a STN temporal constraints are 
used together with disequations in the form x — y ^ r, where r is a real constant, then the problem 
of deciding consistency is still tractable. This extension does not allow the specification of alternative 
constraints but it is interesting because it allows to exclude some solutions maintaining the consistency 
problem tractable. 
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9 Conclusions and Future Work 


In the literature, there are different frameworks and approaches aimed to extend the STN model 
allowing the representation of disjunctive temporal constraints [18,44], but at cost of an exponential¬ 
time consistency check procedure. The only extension with a polynomial time consistency check 
procedure we are aware of is the one of Kumar [40] mentioned in Section 8. 

In this paper, we proposed a novel extension, called Hyper Temporal Network (HyTN), where 
it is possible to represent a new kind of disjunctive constraint, hyper constraint, and to check the 
consistency of a network in pseudo-polynomial time. A hyper constraint is a suitable set of STN 
distance constraints and it is satisfied if at least one distance constraint is satisfied. There could be 
two kinds of hyperarc: multi-head and multi-tail. In a multi-head hyperarc, its distance constraints 
are between a common source timepoint and different destination timepoints. In a multi-tail hyperarc, 
its distance constraints are between different source timepoints and a common destination timepoint. 

A HyTN is said consistent if it is possible to determine an assignment for all its timepoints such 
that all hyperarcs are satisfied. The computational complexity of the consistency problem of a HyTN 
is NP-complete when instances contain both kinds of hyperarc. 

On instances containing either only multi-tail hyperarcs, or only multi-head hyperarcs, the con¬ 
sistency problem can be solved by reducing it, in a very efficient way, to the search of a winning 
strategy in an equivalent Mean Payoff Game (MPG), and exploiting the known winning-strategy 
search algorithms for MPGs. 

Moreover, we presented an empirical analysis of the efficiency of the resulting consistency check 
algorithm. The empirical analysis shows that the proposed algorithm can be effectively used in real 
cases and confirms the general robustness of our approach. 

As future work we are investigating the frontier of practical efficient consistency checking for 
possible generalizations of the HyTN model as, for example, those including contingent constraints [46] 
or conditional ones [45]. 
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